The easy axis of the magnetic response of the conduction electrons of yttrium and 

rare earths 
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We describe a scheme for first-principles calculations of the static, paramagnetic, spin susceptibil- 
ity of metals with relativistic effects such as spin-orbit coupling included. This gives the direction 
for an applied modulated magnetic field to maximise its response. This easy axis depends on the 
modulation's wave-vector q. For h.c.p. yttrium we find the peak response at a q = (0, 0, 0.57)n/c, 
coincident with a Fermi surface nesting vector, to have an easy axis perpendicular to q. This 
is consistent with the helical anti-ferromagnetic order found in many dilute rare earth- Y alloys. 
Conversely, the easy axis for the response to a uniform magnetic field lies along the c-axis. The 
conduction electrons' role in the canting of magnetic moments in Gd — Y alloys and other rare earth 
materials is mooted. 

PACS numbers: 75.40Cx,71.15Rf,75.25+z,71.20.Be,71.20.Eh 
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Conduction electrons play an important role in the 
magnetism of rare earths. Specifically they determine 
the Ruderman-Kittel-Kasuya-Yosida (RKKY) indirect 
exchange interactions which connect the localised 4f mag- 
netic moments of the materials. pj The magnetic re- 
sponse of the conduction electrons, the wave- vector de- 
pendent susceptibility x(q)i governs the ordering of the 
moments. If the greatest value of x(q) is for q = (0, 0, 0), 
the moments order ferromagnetically. A peak for x(q) at 
some finite q, on the other hand, signifies a more compli- 
cated magnetic structure with a modulation wave- vector 
equal to q. In many of the heavier rare earths, e.g. Tb, 
Dy, Ho and Er the incommensurate anti-ferromagnetic 
(AF) ordered states, modulated along the c-axes of these 
h.c.p. crystalline elements, depend strongly on the Fermi 
surface topology of the conduction electrons. The salient 
part of the Fermi surface (FS) that causes this anisotropy 
is the webbing feature that contains flat parallel sheets 
perpendicular to the c-axis forming a strong nesting ef- 
fect HQ. 

The link between FS nesting and the q-dependence 
of x(q) is well-known and has been explored by elec- 
tronic structure calculations. Q There is however an- 
other anisotropy in the interaction between the rare earth 
(RE) magnetic moments that we suggest is also trace- 
able to x(q). This comes from spin-orbit coupling effects 
on the conduction electrons and affects whether the lo- 
calised moments lie perpendicular to the c-axis and thus 
form helical AF magnetic states or instead have compo- 
nents along the c-axis and order into b-axis modulated 
AF phases. Of course the anisotropic crystal field act- 
ing on the 4f-shell plays an important role in this as- 
pect but the similarity between Gd alloys in which the 
crystal field effects are small and alloys of Dy,Tb, etc 
suggests a generic feature from the conduction electrons. 
This is the issue we address in this letter. We describe 
a 'first-principles' theoretical formalism for the magnetic 



response of paramagnetic metals in which all relativistic 
effects such as spin-orbit coupling are included. For the 
first time the easy axis and its dependence upon wave- 
vector q for the magnetic response can be calculated. In 
effect we elucidate a relativistic RKKY interaction. [j| 

We present specific calculations of the magnetic re- 
sponse of the transition metal yttrium. With the same 
h.c.p. crystal structure and electronic configuration as 
the rare earths, apart from the f-electrons, it provides 
an excellent model for their conduction electrons. Much 
is being learnt about the magnetism of rare earths from 
the varied and complex magnetic structures that they 
form in multilayers with yttrium. |6| These range from 
canted, antiphase domain structure in Gd/Y multilay- 
ers to coherent, incommensurate magnetic helices found 
in Dy/Y and Ho/Y structures. The mediating role of 
the conduction electrons of Y for these rare earth ma; 
nctic interactions is crucial in these tailored systems. 
High quality samples of yttrium are more easily come b 
than the heavy rare earths and recent experiments 
have measured directly the 'webbing' feature of its Fermi 
surface that it has in common with several rare earths. 
Moreover the nature of the magnetic ordering of low con- 
centrations of rare earth dopants in Y can be directly 
related to this study. HE El 

We begin by considering a paramagnetic metal sub- 
jected to a small, external, inhomogeneous magnetic 
field, 5h ext (r), and obtain an expression for the induced 
magnetisation Sm(r). We use relativistic density func- 
tional theory (RDFT) ^| to treat the interacting elec- 
trons of the system and derive an expression via a varia- 
tional linear response approach [T3L Phil Il5| . Although 
there are a number of non-relativistic studies of this 
type here we include relativistic effects and pay par- 
ticular attention to the magnetic anisotropy of the re- 
sponse. From our RDFT starting point we make a Gor- 
don decomposition of the current density |12| and retain 
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the spin-only part of the current, namely the spin mag- 
netisation m(r). This results in a spin-only version of 
RDFT. ^3] m which the self-consistent solution of Kohn- 
Sham-Dirac equations is sought, i.e. 

[ ca.p + (3mc 2 + lV ef f [p, m] - (3a .b ef f [p, m] - e ] 

x G(r,r';e) = LS(r - r') (1) 

which describes the motion of a single electron through 
effective fields and a and (3 are Dirac 4x4 matrices. 
G(r, r';e) is the one electron Greens function and the 
charge p(r) and magnetisation densities m(r) can be 
written in terms of it i.e. 



where the x; are measured relative to the positions of 
atoms centred on a;. The non- interacting susceptibility 
of the static unperturbed system is given by 
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m(r) = -Tr.(3a j def{e, p, T) — G(r, r; e) (2) 

where p is the chemical potential, T the temperature 
and f(e,p,T) the Fermi-Dirac function. These expres- 
sions can be converted into sums over fermionic Matsub- 
ara frequencies uj n = i(2n + I^fceT ^t|- The effec- 
tive potential V e ^ [p, m] consists of the usual combina- 
tion of external potential (from the lattice of nuclei), the 
Hartree potential and functional derivative of exchange- 
correlation energy E xc [p, m] with respect to p whilst the 
effective magnetic field h e ff [p, m] is the sum of any ex- 
ternal magnetic field (from magnetic impurities, for ex- 
ample) and the functional derivative of E xc with respect 
to magnetisation. We use the local density approxima- 
tion (LDA) ^5 for E xc . The leading relativistic effects 
contained in the Kohn-Sham-Dirac hamiltonian of eq . (fT|) 
are the well-known mass- velocity, Darwin and spin-orbit 
coupling effects. 

If a small external field 5h ext is applied along a direc- 
tion n with respect to the crystal axes of a paramagnetic 
system, a small magnetisation, <5m(r), and effective mag- 
netic field, 5h e ff are set up. The effective magnetic field 
is given by 8h e 'f [p(r), m(r)] = 8h ext (r) + I xc (r)5m(r) 
where I xc (r) is the functional derivative of the effec- 
tive exchange and correlation magnetic field (within the 
LDA) with respect to the induced magnetisation density. 
The Green's function satisfying equation can be ex- 
panded in a Dyson equation in terms of the unperturbed 
Green's function, G (r,r';e) of the paramagnetic system 
(8h e ff = 0) and perturbation (3a.8h e ff and the first or- 
der terms enable the magnetic response function to be 
obtained. 

For a general crystal lattice with N s atoms located at 
positions a; (I = 1,..,JV S ) in each unit cell, a lattice Fourier 
transform can be carried out over lattice vectors {R;}. 
This can be written 

X n (x ; ,x;, ,q) = Xo( x i, x r ,q) 

JV, 



(3) 



+ ]T / X*(x*,x?,,q)J xc (x?,) 
I" J 

x X fl (x^,4,q)dx^ 



Xo( x ^ x J' »q) = 

G (x;, x[, , k, p + iu> n ) 

VBZ 

x (3a -h G c (x;, , x; , k + q, p + iui n )(A) 

The integral is over the Brillouin zone with wave vec- 
tors k, q and k + q within the Brillouin zone of volume 
Vbz- The sum is over the fermionic Matsubara frequen- 
cies. The Green's function for the unperturbed, param- 
agnetic system containing the band structure effects is 
obtained via relativistic multiple scattering (Korringa- 
Kohn-Rostoker, KKR) theory [l^. We solve equation 
using a direct method of matrix inversion. The full 
Fourier transform is then generated 

X "(q) = (l/V^T,^ 1 ^ 

I I' 

x j dx* j dx{, e^-x;,) x fl (x;j x >, ) qj g> v) (5) 

where V is the volume of the unit cell. Some aspects of 
the numerical methods used to evaluate equations i|3l5|) 
of this type can be found in ^f|. Note that this expres- 
sion for the non-interacting susceptibility can be shown 
to be formally equivalent to one of the familiar type 

rf(q) - /^E'^^ + ?. A) L>- (1 ^ W) (6) 



3,3' 



e,-/(k + q)- ei (k) 



where j is an electronic band index, £j(k) a single 
electron energy, M is a matrix element and /kj = 
/(ej(k), p, T), the Fermi Dirac function. 

The important feature of the response function (equa- 
tions l-il">l) . is its dependence on the direction of the 
magnetic field, n, which vanishes when relativistic, 
spin-orbit coupling effects are omitted. The direction, 
n=(sin#cos<£>,sin#sin<p,cos#) is defined by polar and az- 
imuth angles (0, if). In an h.c.p system such as Y or Sc, 
if an external magnetic field is applied along the n = 
(0,0,1) direction, i.e. the c-axis, where (6, if) — > (0,0), 
Xo is produced. On the other hand \o ^ s the response 
of the system when the field is applied in the a6-plane, 
n = (1,0,0), and (0,ip) -> (vr/2,0). We obtain an 
anisotropy as the difference in the non-interacting sus- 
ceptibility when an external magnetic field is applied 
in two directions with respect to the crystal axes, i.e. 
(x" — Xo)- This is enhanced by exchange and correlation 
effects (cq|3J). The approach presented here is applicable 
to ordered compounds and elemental metals and can be 
modified to study disordered alloys 0] owing to its KKR 
multiple-scattering framework. In order to gauge the im- 
portance of these relativistic effects with atomic number, 
we compare our calculations of the magnetic response of 
Y with its lighter 3d counterpart Sc. 
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FIG. 1: (a) Cross-section of the 'webbing' Fermi surface in 
the H-L-M-K plane for Y. The nesting vector is indicated 
by the arrow. The special point L lies at the centre, (b) The 
enhanced susceptibility for Y and Sc along the T-A direction 
for T = 100K. 

We use atomic sphere approximation (ASA), effective 
one-electron potentials and charge densities in the calcu- 
lations for Y and Sc with experimental lattice constants 
a = 6.89, c = 10.83 and a = 6.24, c = 9.91 respectively, in 
atomic units, arj. [2(| The details of the electronic struc- 
tures using a fully relativistic KKR method compare well 
with those from full potential calculations. [2l| The Fermi 
surface for Y in the H-L-M-K plane is shown in Fig. 1(a). 
It shows two relatively flat parallel sheets. The nesting 
vector qi nc = 7r/c(0, 0, 0.57) is indicated by the arrow. 
Fig. 1(a) is in very good agreement with previous cal- 
culations and experiments 0, . Dugdale et al. Q 
have recently carried out positron annihilation fermiology 
experiments and measured <\ inc = 7r/c(0, 0, 0.55 ± 0.02) 
for Y, Also, Vinokurova et al. 22], measured qi„ c « 
7r/c(0, 0, 0.58). We find a similar Fermi surface for Sc 
with the same nesting vector q inc — 7r/c(0, fLO.57), also 
in good agreement with earlier calculations [2i| ■ 

Our calculations of the enhanced static susceptibility, 
defined by equation © , show a peak at this same wave- 
vector (\i nc = 7r/c(0, 0, 0.57) for both Y and Sc. This is 



FIG. 2: (a) The anisotropy as a function of 9 for Y and Sc 
at q inc = tt/c(0, 0,0.57) and q ~ (0,0,0) and at T = 100K. 
(b) The temperature dependence of the anisotropy of Y and 
Sc at both the nesting vector q;„ c and q . 

shown in Fig. 1(b), where we probe wave- vectors along 
the c-axis from V to A. (The special point A, for h.c.p. 
crystal structures is (0, 0,7r/c).) The temperature is set 
at 100K for these calculations. 

Fig. 2(a) shows the anisotropy Xo (<1> ^> <P) ~ Xofa) as a 
function of 9 for Y and Sc at (\i nc and q Q ~ (0, 0, 0). As 
expected on symmetry grounds, we find the anisotropy 
to be invariant in ip. This means that the magnetic re- 
sponse is insensitive to the direction along which an ex- 
ternal magnetic field is applied in the ab-plane. As 9 
is increased, an anisotropy is observed which reaches a 
maximum at 9 = it/2 as shown in Fig. 2(a). In fact, we 
observe that the anisotropy takes the rather simple form 

X*(q, 0, <p) - X *(q) = (xS(q) " XS(<0) ■ |q x A| 2 (7) 

for h.c.p. crystal structures [24| (|q x n| = s'm8). It is 
therefore sufficient to show X„(q) — Xo(l)i ^° determine 
the easy axes of the magnetic response. Fig. 2(b) shows 
the weak temperature dependence of the anisotropy of 
Y and Sc at both the nesting vector q inc and also at 
q Q ~ (0,0,0). As expected, the anisotropy is an order 
of magnitude larger for Y than for Sc. This difference is 
a result of spin-orbit coupling being more pronounced in 
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the heavier 4d metal Y than in the 3d Sc. We infer that 
a still greater but similar anisotropy should be evident in 
the magnetic response of the conduction electrons in the 
heavier still RE materials. 

It is apparent from Figs.l and 2 that Y shows its 
strongest response at the wave-vector qi„ c which cor- 
responds to the FS nesting and that the response here 
is strongest to magnetic fields directed in the basal ab- 
plane. Here, then, is an explanation as why Gd impu- 
rities, denoted typically as S ions (zero angular momen- 
tum) and thus free from crystal field effects, order into 
helical incommensurate AF states. |l l| It is also a strong 
factor for the similar magnetic order found in other dilute 
rare-earth alloys such as YTb, Y_Dy and YHo. 

Fig. 1(b) shows that yttrium's responses to ferromag- 
netic (x(q) = (0, 0, 0))) and qi„ c -AF modulations are 
close in strength. Figs. 2(a) and 2(b) demonstrate, how- 
ever, that the easy axis for a FM modulation lies along 
the c-axis and therefore contrary to that of the AF one. 
These features are pertinent to the magnetic structure 
of Gd — Y alloys. Pure gadolinium has a ferromagnetic 
phase with the easy axis along the c-axis so that in the 
paramagnetic phase at higher temperatures x(q) peaks 
at q = (0, 0, 0). Moreover its FS here has been observed 
not to possess a webbing feature [2{|. Some recent ex- 
periments 25] have shown that adding more than around 
30% Y to Gd changes the topology of the Fermi surface 
of the paramagnetic phase and the webbing is observed 
again. This is coincident with the alloys' forming helical 
AF states at lower temperatures. Moreover, in this con- 
centration range, neutron-diffraction [j| |2j| has shown 
that application of a modest uniform magnetic field along 



the c-axis leads to a ferromagnetic alignment of the Gd- 
moments along the c-axis. Once the magnetic field is 
switched off the system reverts to its helical AF state. 
This suggests that the alloys' x(q)' s an d anisotropics are 
similar to those shown in Fig. 1(b) and Fig. 2(a) but with 
the relative peak heights of x(q) at q = (0, 0, 0) and q; nc 
to be finely balanced. 

In conclusion, we have described an ab initio theoreti- 
cal formalism to calculate the relativistic static paramag- 
netic spin susceptibility for metals at finite temperatures. 
Since relativistic effects such as spin-orbit coupling are 
included we can identify the anisotropy or easy axes of 
the magnetic response. We applied this formalism to the 
4d metal Y , Its enhanced susceptibility displays a peak 
at the incommensurate wave- vector q mc = (0, 0, 0.577r/c) 
traceable to a FS nesting feature found in the calculated 
electronic structure and fermiology experiments 0, |22| . 
The explanation of magnetic order in rare-earth- Y" alloys 
and superstructures in terms of the FS topology of Y 
is well established both on theoretical [4J and, more re- 
cently, experimental grounds. 8, 25] Results in this letter 
are fully consistent with these findings but we have also 
shown how the electronic structure of Y influences the 
canting of the rare earth moments with respect to the 
crystal axes. The finding that the incommensurate AF 
order of these moments has a helical structure whilst any 
potential ferromagnetism has an easy-axis parallel to the 
c-axis is suggested as a more general property of the con- 
duction electrons of h.c.p. rare-earth materials and has 
a role in the magnetism of RE/Y" multilayers. 
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